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1 Introduction 

Tensor data analysis has successfully developed in various application fields, which 
is useful to seize multi-factor dependence. An nxnxp tensor is a multi-array datum 
T = (Tijk), where 1 < i, j < n and 1 < k < p. A type nxnxp tensor T is denoted 
by T = (Ai, A 2 ; ■ ■ • ; A p ) where Ai denote n x n matrices. An nxnxp tensor T is 
said to be of rank 1 if there is a vectors a = (a 1; a 2 , a n ), b = (b 1 , b 2 , b n ), and 
c = (ci, c 2 , c p ) such that = OjbjC^ for all k, and the rank of a tensor T is 
defined as the minimum of the integer r such that T can be expressed as the sum 
of r rank-one tensors. The maximal rank of all tensors of type nxnxp are also 
defined in obvious fashion and denoted by maxrank(n, n, p).The rank of a tensor 
describes the complexity of a tensorial datum and the maximal rank describes the 
model complexity of a class of tensors of a given type, and so they are very important 
concepts in both applied and theoretical fields. Therefore, the rank and maximal 
rank determination problems have attracted the interest of many researchers, for 
example, Kruskal [12], ten-Berge [24] . and Common et al, [6j, etc. and now have 
being investigated intensively( for a comprehensive survey, see Kolda et al.|llj). 
Atkinson et al.([T] and (2] ) claimed that maxrank(n, n, 3) < 2n — 1. Here we 
introduce an important class of tensors. 

Definition 1.1 A real tensor T = [A\] A 2 ; A p ) is said to be absolutely nonsin- 
gular if x 1 Ax + x 2 A 2 H — • + x p A p is nonsingular for all (xi,x 2 , ...,x p ) 7^ (0, 0, 0). 

Remark 1.2 We called absolutely nonsingular tensors as exceptional tensors in 
Sakata et al. [2T].p2]. 

Sumi et al. [23] proved the claim of Atkinson et al. over the complex number filed 
C without any assumption and proved it over the real number filed 1R except the 
class of absolutely nonsingular tensors. Thus, for the proof of the claim of Atkinson 
et al. over the real number filed M, it is the first thing to determine all absolutely 
nonsingular tensors. Absolutely nonsingular tensors are characterized by the de- 
terminant polynomial defined below. Searching of absolutely nonsingular tensors 
was pursued in Sakata et al. [2"T],|22j in this direction. As well as searching abso- 
lutely nonsingular tensors, the equivalence among them under the rank-preserving 
transformation which is defined below is also important. Note that such equivalence 
relation has also a relation to the SLOCC equivalence of entangled states in the 
quantum communication (for example, see Chen et al. [5]). 



Definition 1.3 For a n x n x p tensor T = (A\\ A 2 \ A p ), the homogeneous 
polynomial in xi, x p of degree n 

fr(xi, X2, x p ) = det(xiA 1 + x 2 A 2 H + x p A p ) (1.4) 

is called the determinant polynomial of a tensor T. 

Then we have the following important characterization. 

Theorem 1.5 If T = (Ai, A 2 ; , . . . ; A p ) is absolutely nonsingular, its determinant 
polynomial f T ( positive definite homogeneous polynomial or negative 

definite homogeneous polynomial. 

Proof Let T = (Ai, A 2 \ , . . . ; A p ) be absolutely nonsingular, and assume that 
there are two points Xq and X\ such that /t(ceo) > and /t(jei) < 0. The line i 
combining the two points Xq and X\ must pass through the origin 0, since in the 
segment [xo,xi] there must be x' such that frix') = and it must be because 
T is absolutely nonsingular. Let take another point x 2 which is not on the line £. 
Then, the line passing x and x 2 does not pass the origin and so f(x )f(x 2 ) < is 
impossible just by the same reason given in the previous sentence. So, f(x )f(x 2 ) > 
0. Next, consider the line passing xi and x 2 , which also does not pass the origin 
and f(xi)f(x 2 ) < 0. This is also a contradiction. After all, there don't exist points 



Xo and X\ such that f(x ) > and f{xx) < 0. This proves Theorem 1.5 I 

It is well known that tensor rank is invariant by typical matrix transformations, 
say, p— , q—, and r— transformations defined below. So, equivalence relation of two 
tensors means that they have a same rank. Thus, to study equivalence among 
tensors is of some importance for rank determination. 

Definition 1.6 For a n x n x p tensor T = (Ai, A 2 \ ■ ■ ■ ; A p ), the following trans- 
formations 

(1) T = (Ai,A 2 ;--- ;A p ) ->■ T = (P A ± ; P A 2 ; ■ ■ ■ ;PA p ) by an n x n matrix 
P G GL(n), 

(2) T = (Ai,A 2 ;--- ;A p ) ->> V = (AQ; PA 2 Q; ■ ■ ■ ; A p Q) by an n x n matrix 
PeGL(n), 

(3) T = (A 1 ; A 2 - ■ ■ ■ ; A p ) -> V = {R ll A l + R 12 A 2 + R 13 A 3 ; R 21 A, + R 22 A 2 + 
R 23 A 3 ; R 31 Ax + R 32 A 2 + R 33 A 3 ) by apxp matrix P 6 GL{n) 

are called as p—, q—, and r— transformations and denoted by T — Y p T',T -^- q 
T'andT — > r T' respectively. Further, if Ti — > p T 2 , the Ti and T 2 are said to be 
in the p— equivalence, q— and r— equivalence are defined analogously. 

Definition 1.7 Let Ti = (A±; A 2 ; ■ ■ ■ ; A p ) and T 2 = (B\\ B 2 ; ■ ■ ■ ; B p ) be two n x 

n x p tensors. If there is a sequence of {T{\ starting from T\ and ending at T 2 , in 
which Ti and Tj + i are in the relation of p—, or q—, or r— equivalence, then Tq and 
Ti are said to be equivalent. 
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Now we can reduce the equivalence relation into a more simple one by the following 
lemma. 

Lemma 1.8 p— , q— and r —transformations are mutually commutative. 

Proof For simplicity, we prove for p = 3, however, the proof is similar for a general 
p. First we prove the commutativity of p— transformation and r— transformation. 
Let 

T\ — > p T 2 — >> T 3 and T\ — Y r T 2 — > p T 3 

We will show that T 3 = T 3 . Let T\ = (A x ; A 2 ; A 3 ) and P = (py) and R = (r^). 
Then, 

T 2 = (PA 1 ;PA 2 ;PA 3 ) 

and 

T 3 = (ruP^i+naP^+ria^a; r 2l P A x +r 22 P A 2 +r 23 P A 3] r 3X PA x +r 32 PA 2 +r 33 PA 3 ) 
On the other hand 

T 2 = (r u A 1 + r 12 A 2 + r 13 A 3 ; r 2X A x + r 22 A 2 + r 23 A 3 ; r 31 A x + r 32 A 2 + r 33 A 3 ) 

and 

Tg = (r n PA 1 +r 12 PA 2 +r 13 PA 3 ; r 21 P A x +r 22 P A 2 +r 23 P A 3 ; r 3l P 'A 1 +r 32 P 'A 2 +r 33 P A 3 ) 

Thus, T 3 = T3, and this means the commutativity of p- and r-transformations. 
The commutativity of q- and r-transformations are proved similarly, p- and q- 



transformations are obviously commutative. This proves Lemma 1.8 



Note that in this paper we consider three cases of (1) P, Q G GL{n) and R £ GL(p) 
and (2) P,Q £ GL(n) and R £ SL(p). and (3) P,Q £ SL(n) and i? £ SL(p). 
The first is called GL(p)-equivalence or simply equivalence, and the second is called 
5 , L(p)-equivalence in short. The third case is called, in a full term, SL(n) x SL(n) x 
SX(p)-equivalence. Lemma 1.8 implies the following theorem. 

Theorem 1.9 Ti and T 2 are GL(p)- equivalent if and only if there is a set of p- 
transformation, q -transformation and r -transformation such that 

T x ^ p T' ^ q T" ->> T 2 

Thus, the equivalence problem of tensors T x = (A^ :, , , : A^) and T 2 = (A^ :, , , : 
A p 2 ^) is reduced to the problem whether the following system of algebraic equations 
for P, Q and R can have a solution or not. 

A (2) = P({J2 r ^ A f ] }Q^ = W»"P ( L1 °) 

These algebraic equations have too many variables to solve even when the size of 
matrices A { is moderate. So, in this paper, we propose to see the problem through the 
determinant polynomial. Then, though we necessarily have to discard the sufficiency 
part of the problem, however, the problem becomes concise and tractable one by 
the following proposition. 
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Proposition 1.11 If T\ and T2 are GL(p)- equivalent, it holds that there is a con- 
stant cGffi and a p x p nonsingular matrix R G GL(p) such that 



f T2 (x) = cf Tl (xR) (1.12) 

So, we can say that 



Proposition 1.13 For two tensors, if the equation \1.12 ) does not hold for any 



constant cGffi and any matrix R G GL{p), they are not GL(p)- equivalent. 



Though the reduced equation ( 1.12[ ) happens to be solved algebraically in some 



cases. However, it is still hard to solve, in general, a system of algebraic equation 
with too many variables. In fact, we need to decide whether a system of ( n+1 H n + 2 ) 
homogeneous equations with (p 2 + 1) variables of degree n have a solution or not. 
So, in this paper, we avoid to solve the problem algebraically and propose to attack 
the problem from a geometric view point, that is, we propose to test non equivalence 
by checking whether the two surfaces of the determinant polynomials of T\ and T 2 
have a same geometric invariants, or not. Here, multi-linear algebra and differential 
geometry intersect through the widow of determinant polynomials. 
The first aim of this paper is to show theoretically that differential geometric invari- 
ants are useful as testers of non-equivalence among absolutely nonsingular tensors. 
The second aim is to show that we can calculate the values of the invariants with 
enough accuracy. Third, we compare the values of invariants calculated by the lat- 
tice method and by the t-design method. And it is shown that the lattice point 
method gives more stable values than the t-design method. As SX(p)-invariant, we 
consider first the volume enclosed by the constant surface and then we consider the 
affine surface area, and thirdly we consider the LP affine surface area of convex body. 
Affine surface area was studied by Blaschke [I] and extended to L p affine surface 
area by Lutwak [20], (also see Leichtwess [13] )• As for a valuation theory of L p 
affine surface area, see the recent papers by Ludwig [15] and Ludwig and Reitzer 
[IT]. Finally, as a general reference of affine differential geometry, see K. Nomizu 
and T. Sasaki [19]. 

This paper is organized as follows. In Section 2, we show how to parametrize the 
constant surface of a determinant polynomial and in Section 3, we review briefly 
some definitions from differential geometry. In Section 4, we argue rough SL(p)- 
invariants. In Section 5, we deal with SX(p)-invariant. In the first subsection, 
we introduce the valuation theory for the set of convex bodies and in the second 
subsection, we argue a volume of the region enclosed by a constant surface as an 
SX(p)-invariant. In the third subsection, we argue the affine surface as a SL(p)- 
invariant. In Section 6, we consider the generalized affine surface, that is, L p affine 
surface area, especially centro-affine surface as a GL(p)-invariant. In Section 7, we 
review the theory of spherical t-design briefly and give a theorem important for ap- 
proximate calculation of our proposed invariants. In Section 8, we give numerical 
values of the invariants calculated by the lattice method and t-design method. It is 
shown numerically that the proposed invariants is usefull to discriminate non equiv- 
alence. In Section 9, the conclusion is given. Finally note that in the following we 
consider mainly the case of n = 4 and p = 3, though some statements are given for 
general n and p. One reason is that absolutely nonsingular tensors are not so easy 
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to obtain for general cases and the second reason is because it is easy to see that 
our method is also available for general cases. The study of much higher values of 
n and p will be given in the future work. 



2 Parametrization of constant surface 

The determinant polynomial of a 4 x 4 x 3 tensor T = (Ai'A^As), ifr{x, y, z) = 
det(xAi + yA 2 + zA 3 ), is a homogeneous polynomial of three variables with degree 
4. We are concerned with the integral invariants of the constant surface dQx — 
{(x,y, z)\fT(x,y, z) = 1} for the special linear group SL(3) and the general linear 
group GL(3). To get such invariants, we need to parametrize this surface by the 
usual spherical coordinate, 

x = r sins cos t = r$ x .(s,t) (2-1) 
y = r sins sin t = r<fr y (s,t) (2-2) 
z = r cos s = r§ z (s, t), (2.3) 

where 0<s<7r,0<t< 2ir. Let x denote the point (x, y, z) on the surface. Putting 
these into the equation f T (x,y,z) = 1, we have 

= — ~r~ — Tj (2.4) 

where 

p(s, t) = f T ($ x (s, t), %(s, t), $ 2 (s, t)). (2.5) 

And so, 

x = P {s\yi* t)j $y(s ' t} ' $z(s ' t]) • (2,6) 



This equation (2.6) gives a parametric representation of the constant surface 8Qt- 



Then, the following is a starting point of this research of the constant surface. 

Theorem 2.7 The constant surface of the determinant polynomial of an absolutely 
nonsingular tensor is a compact set in M. 3 without self-intersection. 

Proof Without loss of generality, we assume that fr{x) is positive definite. If 
x G 8Qt, for any < r < 1 and r > 1, rx is not in dVtx- That is, the constant surface 
is of a star-shaped. This proves that the surface has not any self intersection. Since 
p(s, t) is continuous on the unit sphere it takes a positive minimum and a positive 



maximum. So, x in the equation 2.6 is bounded, which implies the compactness of 



the constant surface. This completes the proof of Theorem |2.7 



The following 8 figures are examples of the constant surfaces of 4 x 4 x 3 absolutely 
nonsingular tensors. 

Note that the numbering of tensors is based on the list of nonsingular tensors 
with elements consisting only of —1,0, 1 found by us. Each figure corresponds to 
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Figure 1: Constant surfaces of the determinant polynomials of 
tensorsFTVol, 3, 10, 20, 99, 119, 207, 237. 

the following tensor and its determinant function respectively. 
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f Tl (x, y, z) = x 4 + 6y 4 + 2z 4 + -3x 3 y - 8xy 3 - 3xz 3 + 5z 3 y + 
7x 2 y 2 + 3x 2 z 2 + 8z 2 y 2 — 2xy 2 z — 8xyz 2 . 
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f n (x,y,z) = x 4 + 3y 4 + 6z 4 -3x 3 z + 2xy 3 + 4y 3 z 

—7xz 3 — 6z 3 y + x 2 y 2 + 5x 2 z 2 — 5z 2 y 2 + 4xy 2 z — x 2 yz + 2xyz 2 . 
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f Tl0 (x, y, z) = x 4 + y 4 + 2z 4 - x 3 y + 2x 3 z + xy 3 + 2xz 3 
—z 3 y — x 2 y 2 + 4x 2 z 2 — 2xy 2 z — 2x 2 yz — xyz 2 . 
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/x 20 (x, y, z) = x 4 + y 4 + 2z 4 - x 3 y + x 3 z + y 3 z - xz 3 + 2z 3 y + 
—x 2 z 2 + 6z 2 y 2 + xy 2 z + x 2 yz + 2xyz 2 . 

1 1 \ / — 1 —1 —1 \ / 1 \ 
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o o i o : 1-110 10 1—1 
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/t 99 (z, 1/, «) = x 4 + 2y 4 + 2z 4 + 3x 3 y + x 3 z + 3xy 3 + 
y 3 z — 3z 3 y + 6x 2 y 2 + z 2 y 2 + Sxy 2 z + x 2 ?/z — hxyz 2 
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f Tll9 (x,y,z) = 
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2x 2 y 2 — 3x 2 z 2 — z 2 y 2 + xy 2 z + x 2 yz + xyz 2 . 
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fT 207 (x, y, z) = x 4 + 2y 4 + 2z 4 - x 3 y - 3x 3 z + xy 3 + 6y 3 z 

— 3xz 3 + 2z 3 y — x 2 y 2 + Ax 2 z 2 + 6z 2 y 2 + 8xy 2 z — 3x 2 yz + 7xyz 2 . 
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3 Notations from differential geometry 

For the use in the following sections, we review some basic notations of differential 
geometry. For more details, for example, see the Nomizu and Sasaki [19] . When 
we denote the parametrized point on the surface by x(s, t), we denote its partial 
derivatives by 

x,(s,t) = (3.1) 
«,(.,«) = ^ 0.2) 
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Definition 3.3 

E=(x 8 ,x 8 ), F=(x 8 ,x t ), G=(x u x t ), (3.4) 

are called the first fundamental coefficients. Putting dx = x s (s, t)ds + x t (s, t)dt, the 
form 

I = (dx, dx) = Es 2 + 2Mdsdt + Ndt 2 (3.5) 
is called the first fundamental form of the surface. 

Definition 3.6 

X S X X t /o _s 

n = - (3.7) 

I \x s x x%\ I 

is called the unit normal vector at the point x(s,t). 
Definition 3.8 Putting 

_ d 2 x(s,t) d 2 x(s,t) d 2 x{s,t) 

Xss ds 2 ' Xss " dst ' Xss ' dp ' 1 

the scalar functions 

L = {x ss ,n) M = (x sU n) N = (x tt ,n) (3.10) 
are called the second fundamental coefficients. The form 

II = -{dX, dn) = Lds 2 + 2Mdsdt + Ndt 2 (3.11) 
is called the second fundamental form of the surface. 

Definition 3.12 At the point P on the surface, let k-y and k 2 be the maximum and 
minimum of curvatures of curves generated by the intersection of the surface with 
the plane spanned by the normal vector and a tangent vector, H — {k\ + &2)/2 is 
called the mean curvature and K = k\k,2 is called the Gaussian curvature. These 
are calculated by 

EN — 2FM + GL LN - M 2 

H = 2(EG-F*) ^ K = EG^- (3 - 13) 

4 Rough 5L(3) invariants 

For checking GL(3)-equivalence between two tensors Ti and T 2 , we need to test the 
equation 

f T2 ( x ) = c f Tl (xR),ce K and R e GL(3). 
Further, S = R/\R\ 1 ^ e SL{3), 

cf Tl (xR) = c\R\ 4/3 f Tl (xR/\R\ 1/3 ) = c'f Tl (xS) with d Gi and S G SX(3). (4.1) 

Thus, GL(3)-equivalence reduces to S'L(3)-equivalence. Then, the following theorem 
holds. 
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Theorem 4.2 For two tensors Ti andT 2 , assume that fr 2 ( x ) — cfr^xS) does not 
hold for any choice o/cGl and anyS* G SL(3). Then, Ti and T 2 are noi GL(3) 
equivalent. 

This justifies to study SX(3)-equivalence among absolutely nonsingular tensors for 
investigating GL(3) equivalence. 

Remark 4.3 If c is negative, then by consider T 3 = PT 2 with \P\ < 0, we have 
f Tz (x) = c f Tl (xR), ceR + and Re GL(3). 

where T 3 is equivalent to T 2 . Thus, we can assume by writing T 3 as T 2 again 
f T2 (x) = c f Tl (xR), ceR + and Re GL(3). 

The following rough SX(3) invariants are useful. 

Theorem 4.4 A convex surface is transformed into a convex surface by a SL(3) 
linear transformation and so, a tensor with a determinant polynomial whose constant 
surface is convex is not equivalent to a tensor with a determinant polynomial whose 
constant surface is not convex. 

Only the tensor of No.l has the convex surface among 8 figures in the Figure 1.1 
and so the tensor of No. 1 is not SL(3) equivalent to all other tensors in the Figure 
1. 

Definition 4.5 A point on the surface is called a singular point if the normal vector 
at the point can not be defined. 

Theorem 4.6 If the constant surface of a tensor Ti has a singular point and the 
constant surface of a tensor T 2 has no singular point, they are not SL(3) (GL{3)) 
equivalent. 

Example 4.7 Let 

/ 1 \ / 10\/0 1\ 
-1000 001 00-10 

1_ 000 -1' -1 000' 0100 
\0 01 / \ -100/ \-1 0/ 

and 

I 1 \ / -1 \ /0010\ 
-1000 0010 0100 

2 ~ 0010 ' 1000 ' 0001 ' 

\o ooi/ \oio o/ yiooo/ 

and the determinant polynomials of Ti and T 2 are given below respectively. Then, 
we have 

f Tl (x,y,z) = (x 2 + y 2 + z 2 ) 2 , 

and 

f T2 (x,y,z) = (x 2 + y 2 ) 2 + z 4 . 
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Note that both of them are positive definite, that is, both of T\ and T 2 are absolutely 
nonsingular. It is clear that the constant surface of fx^x) is a sphere and it has no 
singular point, and on the other hand, that the constant surface of fT 2 { x ) ^ s a conic 
and it has a singular point. Hence, T\ and T 2 are not equivalent. 

Definition 4.8 When we consider a mesh of the parameter space, it produces a lat- 
tice of points on the constant surface. Let K + be the number of lattice points at 
which the Gaussian curvature is positive, and K_ and K be defined in the same 
way. 

Then we have 

Theorem 4.9 The triplet (K + , K_, K ) is an SL (3) -invariant. 

5 SL(3) integral invariants 

The following Figure 2 and 3 shows the figures of convex bodies that are enclosed by 
constant surfaces. The number of figures corresponds to that in our list of absolutely 
nonsingular tensors( Maehara [18J). In this section, we want to find some SL(3)- 



Figure 2: The constant surfaces for Nol, 19, 22, 23, 42, 60, 61, 65 absolutely nonsin- 
gular tensors 

invariant for such convex bodies. For this purpose, the following valuation theory is 
a quite useful. Here, we make a brief summary of the valuation theory from Ludwig 
[15] and [17] . The definition is stated for a general p. 

Definition 5.1 Let K, denote the set of all convex bodies in MP. A functional t(-) 
from K. to M. is called a valuation if it satisfies 



Next theorem is a starting point of characterization of invariant valuation. 




t(K) +t(L) =t(KUL) +t(KHL), K,LeJC. 
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Figure 3: The constant surfaces for No72, 74, 76, 83, 85, 87, 95, 103 absolutely non- 
singular tensors 

Theorem 5.3 (Hadwigerf^j). A continuous valuation t(-) from /C to R is invariant 
with respect to rigid motion if and only if there are constants Cq, C\, c p such that 

t(K) = c Vo(K) + aV^K) + ■■■ + c p V p (K) } (5.4) 

where Vq(K), Vi(K), ■ ■ ■ , and V P (K) are the intrinsic volumes of K. We remark that 
the volumes Vk are called quermassintegrals of K in f73j / and that Vq is the Euler 
index, V p -\ the affine volume of OK , and V p is the volume of the convex body K. 
In the following, we simply denote by a(K) as the affine volume and call it affine 
surface area following [T7\[ , and denoe by V(K) the volume V P (K). 

Definition 5.5 A functional t(-) on IC is said to be equi-affine invariant if it is 
SL(p) -invariant and location invariant. 

The following is essential for us. 

Theorem 5.6 (Ludwig [Lflj ./7^/. /73]/ and JT?}/ ). An upper semi-continuous val- 
uation t(K) from IC to E is equi-affine invariant if and only if there are constants 
Co, Ci, £ R and c<i > such that 

t(K) = c V (K) + Cl V p (K) + c 2 a(K), (5.7) 

where a(K) denotes the affine surface area. 

In short, SL(p) invariant valuation is only the weighted sum of the Euler index and 
the volume V^, (if) and the affine surface area a(K). This means that 

Proposition 5.8 V(K) and a(K) are SL -invariants. 

So, we adopt the volume V(K) and the affine surface area a(K) as indexes of SL(p)- 
equivalence. Further, the next proposition by Lutwak[20j is very useful for us. 

Proposition 5.9 When p = 3,a(K)is homogeneous of degree 3/2, that is, 

a(dK) = d 3/2 a(K), K £ /C . (5.10) 
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5.1 Volume as an 5X(3)-invariant 

We are considering the equivalence relation among absolutely nonsingular tensors. 
As is shown in Theorem 2.7, for such kind of tensors, the constant surfaces of them 
are compact. Note that from Propisition 5.8 the volume of the region enclosed by 
the constant surface is SX(3)-invariant. Then, by the following Gauss's theorem, 
we can calculate the volume by the parametric representation given by the equation 

Theorem 5.11 (Gaussian formula) 

For the region Q enclosed by the space surface dQ, letting fdyAdz+gdzAdx+hdxAdy 
be the differential form of 2nd degree, it holds 

f f fdf do dh\ 

/ fdy A dz + gdz A dx + hdx A dy = / I + + — ) dx A dy A dz (5.12) 
Jan Jn \dx dy dz J 

We denote by V(Q) the volume of the region Q. By this formula, we have 

V(tt) = LHS of the equation (pU2|). (5.13) 

For the present case, by use of the spherical coordinates (s,t), the point of the 
boundary dfl is parametrized as x = r(s,t)($ x (s,t),$ y (s,t),Q z (s,t)). Hence, we 
have 

dy dv 

dy = JLds + -£-dt (5.14) 
ds dt 

:-i/4>^»,+^)«. + 



(-i/4)^?*„ + %^H< 



p 5/4 V ^1/4 



, dz . dz . 
dz = —ds + — dt, 
ds dt 



dp/ds ^ d§ z /ds . , 
dp/dt^ ^ & z /dt 

p5/4 z p 



' 1 /4)-^ 7r ^ + -^)^. 



(5.15) 



Therefore, 



*A* = |(-l/4)^* s + ^)((-l/4)^+^) ds A* 
Similarly, 

*Ad» = ((-l / 4)^. + ^)((-l/4)^. + ^) d .A* 

-1/4)^+*^) ((-l/4)^ + ^) ds A dt , 
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and 

"A* = ((-V4)^ + ^)((-l/4)^* s + ^) rfs A di 

-V4)^*, + ((-V4)^*„ + ^) «. A *. 

By using these, we can calculate the volume of the region enclosed by the constant 
surface of the determinant polynomial. Let T\ and T 2 be two nxnx3 tensors and 
let Qi denote the regions { 

bmx\fTi( x ) < 1}) which are enclosed by the surfaces of {x\f Tl (x) = 1}. Then, by 
SL(3) invariance of volumes, we have 

Theorem 5.16 If Vi ^ V 2 , fr 2 ( x ) fnixR) for any R 6 SL(3), namely, T\ and 
T2 are not SL(n) x SL(n) x SL(3) equivalent. 

For GL invariance, the next lemma is helpful. 

Lemma 5.17 For a determinant polynomial f(x), let V(c) be the volume of Q c = 
{aj|c/(aj) < 1}. Then V(c) = c" 3 / 4 \/(l) for 4 x 4 x 3 case. 

Proof By changing a polynomial f(x) into its constant multiple cf(x), the coor- 
dinates x(s, t) on the constant surface are subject to changes to (^) l l A x(s, t). Hence, 
the integral 

- / xdy A dz + ydz A dx + zdx A dy (5.18) 
3 J an 



is multiplied by c 3 / 4 . This proves the assertion of Lemma 5.17 



Theorem 5.19 Assume thatT\ andT2 be SL(3) equivalent and therefore that there 
is a relation between their determinant polynomials, 



f T2 (x) = cf Tl (xR), 



(5.20) 



where c E R and R G SL(3). Let V\(c) and V^c) denote the volumes of Vtc 
{xlcfxiix) < 1} and Q c = {x\cfT 2 (x) < 1} respectively. Then, it holds that 



(vyy 2 ) 4/3 . 



(5.21) 



Proof The proof is trivial from Lemma 5.17 and omitted. 



From Theorem 5.19 we can know the constant c in the equation (5.20). In the 
next section, it will be made clear that this expresson is helpful for establishing 
GL(3)-equivalence. 
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5.2 Affine surface area as an S'L(3)-invariant 

In this section, for testing SX(3)-equivalence, we propose to use the affine surface 
area, which is an SX(3)-invariant by Theorem 5.8. When p = 3, the affine surface 
area has the following integral expression. 

Definition 5.22 For a smooth convex body K C M 3 , the affine surface area is given 
by 



a(K) 



/ k(K,x) 1 WEG - F 2 dsdt, (5.23) 

JdK 

where n(dK, x) is the Gaussian curvature and E, F and G denote the first funda- 
mental coefficients. 

Next, we show that the affine surface area is useful even as a tester of GL(3)- 
equivalence. Assume that we know the constant c in the relation fT 2 { x ) = cfx^xR) 
with c G M + and R G SL(3) by Theorem 5.19 Then, 

Q 2 = {x\f T2 (x) < 1} (5.24) 
= {x\cf Tl (x) < 1} 
= {x\f Tl (c^x) < 1} 
= c-V'{x\M(x) < 1} 
= c- 1 '^. 

(5.25) 



From Proposition 5^, we have 

a (Q 2 ) = c" 3/8 a(fii). (5.26) 

Thus, we have the following. 

Theorem 5.27 Let T\ and T 2 be absolutely nonsingular tensors. Noting Remark 



4-3, by Theorem 5.19, we can obtain the estimate of c G R + under the assumption 



that their determinant polynomials have the relation fT 2 ( x ) — c /ti(*-R) for some 
unknown constant c G 1R + and an unknown matrix R G SL(3). Then, if a(Q 2 ) 7^ 
c~ 1 / 8 a(Qi) , Tx and T 2 are not GL(3) equivalent. 



By using Theorem 5.17 , this is rephrased as 

Theorem 5.28 Let T\ and T 2 be absolutely nonsingular tensors. Then, if 



a(n 2 ) ^ [y^J G ( fi i)' ( 5 - 29 ) 

T\ and T 2 are not GL (3) -equivalent, where V(Qi) and ^(^2) denote the volume of 
Qi and Q 2 respectively. 
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6 Integral GL(3)-invariant 



In the latter half of the previous section, we presented a procedure to test a non- 
GL(3)-equivalence, however, it is somewhat indirect because we need to estimate 
the constant c before starting the procedure. In this section, we consider a direct 
method handling non-equivalence by using a generalized affine surface area. That 
is, we consider the L q affine surface area, which is an extension of the affine surface 
area and developed by Letwak [20]. Hug [10] gave an equivalent definition. The 
following is the Hug' s definition. 

Definition 6.1 

L q a{K) = / K (K,x)^ida K (x) (6.2) 

JdK 

where 

k(K,x) 

k (K,x) = -. ; ' (, — . 6.3 

and d(TK{x) is called a cone measure defined by 

daxix) = (x,n(K,x))dx, (6.4) 
and n(K, x) denotes the outer normal at x on dK. 

When q = 1, L q (K) becomes the affine surface area a(K), and when q = p, it 
becomes a classical centro-affine surface (K) thta is defined as 

a c (K)= [ K (K,xy/ 2 da K (x), (6.5) 

JdK 

which is known to be GL(p)-invariant. The characterization of a general GL(p)- 
invariant functional is given below. 

Theorem 6.6 (Tudwig and Reitzsner [T7] ) Let )Cq be the space of convex bodies 
that contain the origin in their interiors. An upper semi- continuous functional £(•) 
from /Co to 1R 1 is GL(p) -invariant if and only if there are nonengative constants Co 
and C\ such that 

t(K) = c V (K) + Cl a c (K). (6.7) 

7 Spherical design 

According to our experiments, the numerical integrations of the invariants must 
be accurate at least 2 decimals. So, the calculations of the invariants are a little 
bit heavy. In this section, we consider the t-design method as an substitute of the 
nemerical integrations. The spherical design was initiated by Delsarte et al. [7] and 
has been studied by several researchers, for example, see Bannai and Bannai [3]. It 
is defined as follows. 
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7.1 An overview of spherical design 

Definition 7.1 A finite set X on the sphere is called t-spherical design if the fol- 
lowing equality holds that for any polynomial f(x,y,z) with a degree less than or 
equal to t, 



f(x, y, z)da = ^ f{x,y 



(7.2) 



(x,y,z)ex 



where S 2 denotes the unit sphere of M 3 and dadenotes the surface element of the 
sphere and \S 2 \ denotes the surface area of the sphere. 



A parametrized integral formula of the equation |7.2| is given by 

JV 



/ f(s,t) sm(s)dsdt = ^-^2f{si,ti) 

i=i 



(7.3) 



where (si,ti),i = 1,2, ....,N are the corresponding parameters to the design points 
in X. One point to overcome for our purpose is that we need to integrate some 
nonlinear functions that are not polynomials and hence we can not use any t-design 
directly. However, we can rely on the next theorem to solve this point. 

Theorem 7.4 Let f(x,y,z) be an continuous function over the unit sphere and let 
€t be a positive number such that 

\f(x,y,z) -p(x,y,z)\ < e t 

uniformly for some polynomial p(x,y,z) with degree less than or equal to t. Then, 
it holds that 



1 N 

/ f(x,y,z)dS - An— y2f(xi)\ < 8ne t 
JdS ^ i=1 



(7.5) 



Proof 



< 



/ f(x,y,z)dS - 4n^-^2f(x l ) 

JdS i=1 
f 1 N 

/ f{x,y,z)dS-4n— 2Zp{xi) 
Jas ^ i=1 



(7.6) 



l N x TV 



i=l 



1=1 



f(x,y,z)dS - / p(xi)dS 



as 



| ds 



+ 



Air 
TV 



N 



i=l 



< J e t dS + 4ire t 

= 8TT6 t 
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Remark 7.7 By the above theorem, we need not to know the best approximate poly- 
nomial concretely in order to obtain an approximate value of the integration, and it 
is enough to use f(x,y,z) itself. Moreover the error of the approximation is bounded 
from above by the multiple of e t by 8rr. For a substantial evaluation of the approx- 
imation, we need to know e t . The problem is interesting, however, it is a little bit 
heavy task at present, and so it is postponed to the future work. 



7.2 Calculation of integral invariants by a 20-design 

Using the result of the previous subsection, we consider the integration 

a[K) = J K(s,t) 1/4 VEG- F 2 dsdt. (7.8) 

where s, t moves < s < ir,0 < t < 2n. This integration can be thought to be an 
integration over the unit sphere by 

a(K) = [ K(s,t)^WEG - F 2 dsdt (7.9) 



y/EG - F 2 

I K( S ,t)^ VEG ~ F -dS, 
JdS sm s 



= / /t^t) 1 / 4 ; sin sdsdt 

J(s,t,)e[0,n]x[0,2ir} S m S 



(7.10) 



where dS = sin(s). Hence, 



p(x,y,z) = K(8,t) 1 ' 4 VEG F " (7.11) 

sm s 

is taken to be a function over the unit sphere and so the integral invariant can be 
approximated by the right hand side of the equation below. 

/ k(s, t) 1/4 VE ° ~ F -dS ~^jr p( Xi , Vi , Zi) (7.12) 
JdS sms iV i=1 



The values of invariants calculated by the lattice method and the 20-design method 
will be give in the next section. The 20-design method show very nice approxi- 
mations in some cases, however, do not show good approximations for other cases. 
That is, for our integration of invariants, the spherical design method does not give 
stable values, unfortunately. This might suggest that we need to use design with 
more higher degree than 20. 



8 Effectiveness of the invariants as testers of non- 
equivalence 

In this section, we will show the effectiveness of the numerical values of the invariants 
as testers of non-equivalence. We numerically calculated the volume V(Q), the affine 
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surface area a(f2)and centro-affine surface area a c (Q) of the region Q = {x\f(x < 1} 
defined by the determinant polynomials fr{x). As examples, we calculate them for 
the 16 tensors which are in /Co, whose constant surfaces are figured in Figures 2 and 
3 in the section 5. The numerical calculations are performed in two way, that is, 
by the lattice method and by the t-design method, and they are compared. As for 
the t-design method, we use the 20-design named des.3.216.20 in [8] which has 216 
points. In the tables below, M1-P2-G5, M6-P2-G, M1-P2-G7 and 20-design denote 
the globally adaptive integration with accuracy of 5 digits, pseudo-Monte Carlo 
integration, the globally adaptive integration with accuracy of 7 digits by 64 decimal 
calculation and 20-design method by IEEE754 decimal calculation, respectively. For 
all calculation were done by Mathematica. Table 1 shows that the SL invariance of 



Tensor 


VO 


VI 


V2 


V3 


T001 


2.9197794095194 


2.9197794099529 


2.9197794089308 


2.9197794061274 


T019 


4.0314824331814 


4.0314824340674 


4.0314824332515 


4.0314824319603 


T022 


3.6306602017309 


3.6306602004447 


3.6306602054741 


3.6306602016552 


T023 


3.4355628950802 


3.4355628819358 


3.4355628878857 


3.4355628897838 


T042 


3.7515624235646 


3.7515624142272 


3.7515624197586 


3.7515624152774 


T060 


2.1440485535226 


2.1440485507771 


2.1440485551454 


2.1440485550215 


T061 


2.8594583429857 


2.8594583441125 


2.8594583445567 


2.8594583445567 


T065 


3.1084258968340 


3.1084258946417 


3.1084258957994 


3.1084258984271 


T072 


4.6861403575076 


4.6861403597489 


4.6861403542060 


4.6861403560079 


T074 


3.6302252513670 


3.6302253269919 


3.63022533280632 


3.6302253350968 



Table 1: Volumes by M1-P2-G7 :Tn, where n = 001,0019,022,023,042,060,061, 
065, 072 and 074 Each line denoted as Tn-0 lists the value of the original tensor and 
the lines Tn-i,i = 1,2, ...,5 list the values for the transformed tensors of Tn by a 
randomly chosen matrix of SL(3). 

volumes of the redions enclosed by the constant surface is clealry seen numerically 
for every absolutely nonsingular chosen tensors. Tables 2 and 3 of the affine surface 
area show that the affine suraface area is SL invariant and that all relevant tensors 



are not SX(4) x 5X(4) x SL(3) equivalent mutually. From Theorem 5.28, combining 
the volume data, we also conclude that they are not GL equivalent. This last fact 
is also derived by a direct usage of the centro-affine surface data which is seen in 
Table 4 and Table 5. 

Indeed, Tables 4 and 5 show that the centro-affine surface area is really GL 
invarinat, and that three point decimal accuracy will be sufficient to detect non 
GL(3)-equivalence between 4x4x3 absolutely nonsingular tensors, whose elements 
consists of only -1,0,1. The M1-P2-G7 method seems clearly the best for discrimi- 
nating the tensors relating to GL nonequivalence. 



18 



1 on cnr 


mi P9 


Mfi P9 C ^ 


M1 P9 C7 

1V1 1-JT Z _ V_T f 


Ofl nDcirrn 
ZrU-Llcblg,!! 


rnnn i r\ 

1UU1-U 


9.961493457 


9.962 /96404 


a ap 1 /i T1 /Iao 

9.961471493 


r\ nno 1 T 

9.90317 


1001-1 


9. 961470358 


9.961249135 


n n/^ 1 a ^7 1 /i on 

9.961471489 


o 7onno 

8.73023 


Tnm 9 

1 UU1-Z 


Q Qfi1 A 71 1 Q 1 ! 

y.yoi4 i noo 


Q Q71 9fifi7^n 

y .y ( izoo i ou 


Q Qfi1 /I 71 A 8fi 
y.y014 ( 14oO 


y.yoozo 


i uu i-o 


Q Qfi1 47n f ?97 
y.yui4 i uoz i 


q Q c iQ c ;nQ7nQ 
y .yoyouy i uy 


Q Qfi1 471 474 
y.yui4 1 14 < 4 


Q 7Qn^7 
y . i yuo i 


Tnm a 

1 UU l -z ± 


Q Qfi1 471 1 8fi 
y.yui4 i iiou 


Q Q7Q4^fi1 8fi 
y.y / y4ouiou 


Q Qfi1 471 478 
y.yui4 1 14 / o 


1 n 88^fi7 

lU.OOOU ( 


Tnm ^ 
1 uu i-o 


Q Qfi147499n 
y.yui4 / 4zzu 


Q QQ71 SnQSQ 
y .yy / louyoy 


Q Qfi14714Qn 
y.yui4< i4yu 


1 n QQ978 
lu.yyz i o 


1019-0 


11.560007113 


11.560055546 


11.560007991 


1 1 70^77 

11.87277 


rpn 1 r\ 1 

1019-1 


11. 560007742 


1 1 rrnm ^7ono 

11.552017302 


1 1 r/'nnn^nno 

11.560007993 


1 1 /"'noon 

11.69239 


Tm Q 9 
i uiy-z 


11 ^nnnns^s 

ll.OOUUUoOoo 


1 1 ^Q8fifi/I9/1 

i i.ooyoou4Z4 


1 1 ^Rnnn7QQ'? 
1 1 .oouuu i yyo 


11 c;i *5 /i /i 

1 1.0 1044 


Tm Q 1 


1 1 ^finnn7fiQ9 

ii.ouuuu i uyz 


1 1 ^m nnQQ7i 

ii.ouiuuyyi i 


1 1 ^nnnn7Q8Q 

ii.ouuuu i yoy 


io.ou i uy 


Tm Q 4 


1 1 ^RnnnszLQzL 


1 1 ^4^9nnm 7 


1 1 ^nnnn7QQi 

ii.ouuuu i yyi 


1 n 4n9n^ 


Tm q ^ 


1 1 ^nnnm q9zl 
ii.ouuuuiyz4 


1 1 ^81 7fi^99 
1 l.OOOl 1 UOZZ 


1 1 ^Rnnn7QQR 
1 1 .ouuuu i yyu 


1 1 4Q^Q^ 
i i.4yoyo 


1022-0 


11.020675684 


11.024551464 


1 1 AOA/^^7/1 1 OC 

11.020674135 


1 1 A CO AO 

11.05202 


1022-1 


1 1 non/"7000i 

11.020673831 


11 ni /^n 1 ^7 ,1 o A 

11.016947424 


11.020674138 


11 A CT O A CT 

11.45345 


1 \JZZ-Z 


1 1 n9nfi7fii oc; 
ll.UzUtw olyo 


1 1 n97c;9c;Qc;n 
11. Uz 1 ozoooU 


1 1 n9nfi7/i 1/17 


1 l.o4ooo 


Tn99 1 


1 1 n9nfi7 f ?91 4 

ll.UZUU i OZ14 


1 1 m nnnnn^Q 
i i.uiuuuuyoy 


1 1 n9nfi74i 4n 

1 1 .UZUU ( ^il^iU 


1 1 n7^94 

1 l.U i OZ4 


Tn99 ZL 
1 UZZ-4 


1 1 n9nfi7^Qci 
n.uzuu ( ooyy 


1 1 n^l ^QfiQ^ 
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1023-0 


1 A ^Ti a /i on 

10.771760482 


10.773095422 


10.771760351 


10.73293 




10.771758801 


i n '7'7 / i '7crr\0/^ C 

10.774759865 


10.771760349 


9.26881 




1 n 771 7^Q^m 
J-U. i i 1 1 oc/oui 


1 n 79^&4^9Q1 


1 n 771 7nn^9 

1U. I 1 1 1 uuooz, 


1 n Q4^87 

1 U. VHiOO t 


T023-3 


10.771759730 


10.766135806 


10.771760352 


13.23848 


T023-4 


10.771757516 


10.773059224 


10.771760350 


10.78732 


T023-5 


10.771759533 


10.773749461 


10.771760351 


10.88149 


T042-0 


11.136697741 


11.136755128 


11.136697332 


10.99424 


T042-1 


11.136695637 


11.140565725 


11.136697314 


11.28015 


T042-2 


11.136699257 


11.140835239 


11.136697323 


11.89052 


T042-3 


11.136721676 


11.203272583 


11.136697308 


12.13058 


T042-4 


11.136696147 


11.106329119 


11.136697270 


9.81007 


T042-5 


11.136697731 


11.107102150 


11.136697313 


13.99432 



Table 2: AfRne surface area:Tn, where n = 001, 019, 022, 023 and 042 Each line 
denoted as Tn-0 lists the value of the original tensor and the lines Tn-i,i = 1, 2, 5 
list the values for the transformed tensors of Tn by a randomly chosen matrix of 
SX(3). 
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10.55820 


T074-2 
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10.732331733 
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10.57017 



Table 3: AfRne surface area:Tn, where n = 060, 061, 065, 072 and 074 Each line 
denoted as Tn-0 lists the value of the original tensor and the lines Tn-i,i = 1, 2, 5 
list the values for the transformed tensors of Tn by a randomly chosen matrix of 
SX(3). 
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Table 4: Centro-affine surface area:Tn, where n = 001, 019, 022, 023 and 042 Each 
line denoted as Tn-0 lists the value of the original tensor and the lines Tn-i,i = 
1,2, ...,5 list the values for the transformed tensors of Tn by a randomly chosen 
matrix of GL(3). 
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Table 5: Centro-affine surface area:Tn, where n = 060,061,065,072 and 074 Each 
line denoted as Tn-0 lists the value of the original tensor and the lines Tn-i,i = 
1,2, ...,5 list the values for the transformed tensors of Tn by a randomly chosen 
matrix of GL(3). 
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9 Conclusion 



We treated the SL{4) x SL(4) x SL(3), GL(4) x GL{A) x SL(3) or GL(4) x GL{A) x 
GL(3) non-equivalence problem of 4 x 4 x 3 absolutely nonsingular tensors. We 
proposed a method to addres to the problem through the determinant polynomials. 
Furthermore we proposed to solve the problem by differential geometric SL(3) or 
GL(3) invariant of the constant surface of the determinant polynomials. From the 
numerical analysis by Mathematica, it was shown that the stable values of invariants 
are obtainable numerically and also it was shown that the affine surface area and the 
centro-affme surface area are useful to detect the non-equivalence. This means that 
the algebraic problem: whether a system of algebraic equations with many variables 
can have real solutions or not, can be resolved by differential geometric methods. It 
is a nice link between algebra and differential geometry. Second, we investigated the 
spherical design method for calculating invariants. At present, we think that the 
values given by the adaptive lattice methods are more reliable than those given by 
the spherical design method. In some future work, we expect to extend the result 
to more higher dimensional tensors and to know why the spherical design method 
does not give stable values of invariants. 
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